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Recent experiments indicate that Na^^Os, a material in which s=l/2 Ir local moments form a 
three dimensional network of corner-sharing triangles, may have a quantum spin liquid ground state 
with gapless spin excitations. Using a combination of exact diagonalization, symmetry analysis of 
fermionic mean field ground states and Gutzwiller projected variational wavefunction studies, we 
propose a quantum spin liquid with spinon Fermi surfaces as a favorable candidate for the ground 
state of the Heisenberg model on the hyper-kagome lattice of Nadirs Og. We present a renormalized 
mean field theory of the specific heat of this spin liquid and also discuss possible low temperature 
instabilities of the spinon Fermi surfaces. 



Introduction: Na^^Os is a recently discovered three 
dimensional (3D) frustrated quantum magnet [l[ . The Ir 
atoms in this insulating compound have s=l/2 local mo- 
ments and form a 3D network of corner sharing triangles 
called a 'hyper-kagome' lattice a cubic lattice whose 
unit cell is shown in Fig. [TJ High temperature magnetic 
susceptibility (\) measurements in this material suggest 
that the Ir moments have strong antiferromagnetic cor- 
relations with a Curie- Weiss temperature &w ~ —650^. 
The observation of a large x an( i entropy at low temper- 
ature indicates that gapless spinful excitations survive 
for T <C <dw At the same time, x an d specific heat 
measurements reveal no sign of any magnetic order or 
any other symmetry breaking down to T ~ 0.5K, nearly 
three orders of magnitude lower than Qw, suggesting 
that Na4lr308 may be the first example of a 3D quan- 
tum spin liquid which does not order down to T = 0. It 
joins a small but growing list of recently discovered frus- 
trated s=l/2 quantum magnets Q which appear to have 
quantum disordered ground states possibly supporting 
fractionalized excitations. 

These experiments motivated a study of the classical 
Heisenberg antiferromagnet on the hyper-kagome lattice 
0. This model was found to order into a coplanar 'clas- 
sical nematic' state at low temperatures, T < J/1000, 
where J is the nearest neighbor antiferromagnetic ex- 
change coupling. However, quantum effects are clearly 
significant at such low temperatures. A subsequent study 
of the quantum Heisenberg model, using an Sp(N) mean 
field theory, uncovered a candidate quantum spin liquid 
ground state with topological order as well as an in- 
teresting magnetically ordered state which is proximate 
to this spin liquid [4]. However, this 'bosonic' spin liquid 
has a nonzero spin gap which is at odds with recent ob- 
servations, that gapless spin excitations survive down to 
T ~ 0.5K [5(, unless the spin gap is anomalously small. 
Another difficulty of this proposal is that there should be 
a finite temperature transition from the Z2 spin liquid to 
the higher temperature paramagnetic phase while there 
is no clear signature of a phase transition in thermody- 
namic measurements [l|. 



FIG. 1: Exact diagonalization (ED) results [solid line] for the 
inverse uniform magnetic susceptibility l/\ compared with 
experiments [open circles]. The ED was done on a single unit 
cell of the hyper-kagome lattice (inset) with J = 304K chosen 
to reproduce the high temperature experimental 1/x- 



Here we pursue a completely different line of attack 
and attempt to build a 'fermionic' spin liquid theory of 
the hyper-kagome Heisenberg model. This formulation 
has the virtue that gapless spin liquids emerge as sta- 
ble phases at mean field level and beyond without any 
need for fine tuning 0, @] • The main results of our paper 
are as follows, (i) We find that of a number of candidate 
spin liquid ground states we have explored, a particularly 
simple fermionic spin liquid state, one which supports 
Fermi surfaces of spinons, emerges as a promising candi- 
date for the ground state of the nearest neighbor Heisen- 
berg model on the hyper-kagome lattice. This result is 
obtained by a combination of exact diagonalization, a 
projective symmetry group (PSG) analysis Q of mean 
field ground states, and Gutzwiller projected variational 
wavefunction calculations, (ii) We then show, using a 
Gutzwiller renormalized mean field theory [§], that the 
specific heat of this spin liquid state is quite similar to 
the experimentally observed specific heat of NaglrsOs for 
T > 5K. This spinon Fermi surface state therefore seems 
to be a good starting point to understand the physics of 
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this material over a wide range of temperatures in the 
same way that Fermi liquid theory is a good starting 
point to understand conventional metals. However, as 
in conventional metals, the Fermi surfaces could be un- 
stable, at very low temperature, due to small additional 
interactions, (iii) We considered a symmetry analysis of 
possible low temperature instabilities of the spinon Fermi 
surface state. Some of the states resulting from this anal- 
ysis support line nodes for spinon excitations. We discuss 
implications of such line node states for the specific heat 
data. 

Model and exact diagonalization: We begin with an 
exact diagonalization (ED) study of the nearest neighbor 
s = 1/2 Heisenberg model 



field Hamiltonian 



H — — Jij @i ■ Sj , 



(1) 



on a single 12-site unit cell of the the hyper-kagome 
lattice formed by the Ir sites in Na^^Os (see inset of 
Fig. [T]) . Here Jy is the exchange coupling on the bond 
ij, and we keep only the nearest neighbor antiferro- 
magnetic exchange interaction, J > 0. Fig. [T] displays 
the ED result for x" 1 ^) with a choice of J = 304.fr 
which, as shown, reproduces the experimental data at 
high temperatures. We find that for T = 200-300A', the 
X -1 (T) from ED can be fit by a "Curie- Weiss" law with 
Q w « — 730K . This is only an apparent Curie- Weiss be- 
havior in the sense that it is only valid for a limited range 
of temperature around J and not for T» J where one 
recovers the usual Curie- Weiss law with 6^^ J = — z J/4, 
with the coordination number z = 4 on the hyper-kagome 
lattice. The upturn in x^{T) in the ED for T < 50K 
arises from a nonzero spin gap on a single unit cell. 

Experiments observe a broad peak in C /T at T p « 
25K with a peak height (C/T) max w 55mJ/K 2 /mol-Ir 
While finite size effects are clearly important in the 
ED calculations at low temperatures, it is nevertheless 
encouraging that the ED result for C/T of the Heisenberg 
model with J = 304K shows a broad peak at T p w 2QK 
with a peak height of (C/T) max w 70mJ/if 2 /mol-Ir. 

Since the s = 1/2 Heisenberg model appears to capture 
some aspects of the experimental data on Na^raOg, we 
now turn to an analysis of possible fermionic spin liquid 
candidates for the ground state of this model in order to 
see if we can understand the emergence of a gapless spin 
liquid ground state in Na^^Os- 

Hyper-kagome spin liquid states: We begin by repre- 
senting spin operators in terms of fermionic spinors with 
the constraint of one fermion per site. The Hamiltonian 
H can then be decoupled in both the hopping (Hartree- 
Fock) and pairing (Bogoliubov) channels within a mean 
field approximation. This decoupling leads to the mean 



(ij) 



(2) 



written in a manifestly SU(2) invariant form where ^S>J = 



(fih fi\) T i s a Nambu spinor, 



_ ( Xij 
A*. 



'Xij 



(3) 



are the mean field hopping and pairing amplitudes and Si 
is a lagrange multiplier which on average enforces, in an 
SU(2) invariant manner, the single occupancy constraint 

(fiofia) = 1- 

We seek the ground state of H. Guided by experiment, 
assume the ground state is a spin liquid, that it does 
not break any lattice symmetries, global spin-rotation 
symmetry or time-reversal symmetry. Remarkably, even 
aside from the known SU(2) gauge redundancy, imposing 
all these symmetries on the mean field ground state does 
not single out a unique choice for f/y and a*. The identifi- 
cation of non-redundant spin liquid phases, therefore, re- 
quires careful examinations of the combined space group 
and gauge transformations. These fundamental symme- 
try relations, proposed by Wen 0] and dubbed "projec- 
tive symmetry group (PSG)" can be used to classify all 
possible symmetric spin liquid phases. Under a PSG 
transformation, spinors transform as ^ — ► • ^x(i) 
where Gf is an SU(2) gauge transformation associated 
with the space group transformation X. Transforming a 
mean field ansatz by: 



Ui. 



G 



A' 



i ■ t —> Gf (a X (i) ■ r)Gi 
(4) 

then leaves Hmf invariant. Naturally, these PSGs are 
subgroups of the combined space and gauge groups, 
which leads to strong constraints on the possible choices 
ofGf. 

We have constructed a systematic classification of spin 
liquid ground states by constructing all PSGs with non- 
trivial Gf associated with the point group of the hyper- 
kagome lattice. This group turns out to be equivalent to 
the octahedral group O and consists of two-fold rotations 
about each site, three-fold rotations for each triangle and 
four-fold screw rotations for each thread (see Ref. Q 
for a discussion of threads). Details of our calculation 
will be presented in a future publication [11] . Here, for 
simplicity, we focus on the family of states 



-fJ,Z 



(5) 



where Xij is rea l an d positive, Ay is real but alternates 
sign as discussed below and only the bonds ij that have 
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a finite exchange Jij have finite fTy. This family of states 
covers most of the states resulting from our PSG analysis. 

From a symmetry perspective, Xij an d Ay are cho- 
sen to be invariant under translations and three-fold ro- 
tations through each triangle. However, two-fold rota- 
tions about each site and the four-fold screw rotations 
both need to be followed by the gauge transformation 
Gf — ir3, where X is either of these transformations. 
This second requirement fixes the sign of the pairing 
fields Ajj . In addition to these spatial symmetries, we 
have imposed time reversal (T) invariance by requiring 
that a T transformation followed by Gj commute (or 
anticommute) with the spatial transformations. Since 
T sends f7 y — » —Uij, we found Gf — ir 2 satisfies all 
requirements. The combination of all these symmetries 
completely determines the form of Uij in Eq. [5] 

It turns out that due to enhanced symmetry in spe- 
cial limits, the ansatz of Eq. ([5|) describes three differ- 
ent spin liquid states: the U(l)-uniform state, the U(l)- 
staggered state and the Z2 state. The U(l) uniform state 
has Xij > 0, no pairing (A y - = 0) and a U(l) phase in- 
variance. On the other hand, the U(l) staggered state 
has no hopping Xij = an d finite Ay ^ which alter- 
nates sign on adjacent triangles. In this state, the U(l) 
phase invariance can be understood after noting that an 
SU(2) gauge transformation can rotate this pure pairing 
state into a pure hopping state (with hopping x%j — Ay ) . 
Lastly, the Z2 state arises when both pairing and hopping 
are present. 

In general, we need not keep the time reversal symme- 
try of the ansatz in Eq. |5]). If we let Xij — u ij cos0y 



and Ai 



lij sin 6ij, so that sign( 



;') = sign(Ay), we 



can extend the ansatz to 



Uij = iUij exp{- 



r}, 



-[iZ 



(6) 



where hij — z cos 9ij + x sin 0y . This extended form then 
has all the same spatial symmetries of Eq. ([5]) but recov- 
ers time reversal invariance only at v = 7r. 

Having discussed the fermionic projective symmetry 
group approach to construct a set of distinct spin liquid 
states on the hyper-kagome lattice we would like to know 
which of them are viable candidates for the ground state 
of the Heisenberg model. In order to address this issue, 
we compare the energies of these different states. 
Energetics of candidate spin liquid states: We have 
computed the ground state energy for the above class 
of states in mean field theory as well as by a numeri- 
cal Gutzwiller projection of the mean field states which 
yields a physical spin wavefunction. The Gutzwiller pro- 
jected energy is computed using the variational Monte 
Carlo (VMC) method @. 

• U(l)-uniform state: The mean field ground state en- 
ergy per spin is E^ f — —0.144 J. After Gutzwiller pro- 
jection, we find a variational energy E^°f « —0.424 J, so 
that E^°f/E'^ i{ w 3. The energy of the projected state 
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FIG. 2: Energetics of the Z2 family of states parameterized 
by 9. Points show VMC data and the solid line is the mean 
field result multiplied by gj w 3. The uniform state (9 — 0) 
appears to have the lowest energy in this class of states. 



compares favorably with the ED result on a single unit 
cell, E cd = — 0.454J. In the pre-projected state, three 
spinon bands cross the Fermi level. One Fermi surface 
is electron-like and centered at K = (0,0,0), while the 
other two are hole- like and centered about K = (tt, tt, tt). 
All three have Uf ~ 0.2-7T. 

• U(l)-staggered state: The mean field energy of this 
state is -E^g = — 0.122J. Due to flat bands at the chem- 
ical potential in this state, the energy of the projected 
wavefunction depends somewhat on our selection of the 
subset of the flat band states we fill with fermions in 
the preprojected state. For various choices that we have 
explored the estimated VMC energy is about E^°^ ~ 
— 0.37J, significantly higher than the uniform state. 

• Z2 state: As seen from Fig. [21 the mean field energy 
of the Z2 states parametrized by 9 = 1 9^ | is higher than 
that of the U(l) uniform state (which corresponds to 9 = 
0). Even after projection, the U(l) uniform state appears 
to have the lowest energy, although the energy is quite 
flat as a function of 9 for 9 < 0.1 tt as in the mean field 
theory. The projected energy differs from the mean field 
value by about a factor of three over a wide range of 9. 

• Chiral states: We have also checked the energetics of 
the time-reversal symmetry broken chiral U(l) spin liquid 
ansatz. The uniform U(l) state is stable against such 
symmetry breaking. However, the staggered U(l) state 
energy is lowered by breaking time reversal symmetry 
Nevertheless, the lowest energy obtained in this manner, 

E IkJ 0.39J, is still higher than the uniform U(l) 

state energy. 

In summary, the U(l) uniform state appears to be the 
most favorable candidate for the ground state of the near- 
est neighbor Heisenberg model on the hyper-kagome lat- 
tice. However, as seen from Fig. [2j the energy is a rather 
flat function of 9 for small values of 9 < O.I71-. We there- 
fore cannot rule out the possibility that further neighbor 
couplings or an extended variational ansatz with more 
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variational parameters will not favor such a Z<i state with 
a small pair amplitude. We discuss this further in our 
concluding section. 

Application to the specific heat of Na^rsOs: Mo- 
tivated by our variational ground state calculations we 
next turn to specific heat of the uniform U(l) state in 
order to compare with the data on Na^^Os- Since 
we cannot implement the Gutzwiller projection exactly 
for computing finite temperature properties in any sim- 
ple manner, we will use a renormalized mean field the- 
ory (RMFT) H to make progress. The RMFT analy- 
sis suggests that a part of the effect of projection can 
be taken into account on top of the mean field the- 
ory by Gutzwiller renormalization factors. For instance, 
(Si ■ Sj) pxo j = g.j(Si ■ Sj) m { defines the renormalization 
factor gj for the energy. This renormalizes the mean field 

spinon bandwidth. From our calculations, we find that 
E proi/ E nd K 3 which imp ii es gj ~ 3. 

We next set gj — 3 and compute the quasiparticle con- 
tribution to the specific heat of the U(l) uniform state 
in the RMFT. Fig. [3] shows the heat capacity computed 
this way. We emphasize that this is a zero-parameter fit. 
As seen, the overall behavior of C/T is in broad agree- 
ment with the data. Remarkably, for 5K < T < 25K, 
we find C /T shows a strong, almost linear, T -dependence 
similar to experiment arising simply from the spinon dis- 
persion in the uniform U(l) state. However, as expected, 
the computed C/T eventually saturates below about 5K 
due to the spinon Fermi surfaces, leading to a nonzero 
7 ks 10 mJ/K 2 /molTr. A more precise estimate of 7 
requires the projection of excited states — this is numer- 
ically complicated due to the many bands in this system 
and was not attempted in this work. We emphasize that 
the broad agreement with experimental data already pro- 
vides a nontrivial check of our theory. Integrating C/T 
we also present a comparison between the entropy of this 
state which also shows reasonable agreement with the 
higher temperature data. 

Discussion: We have argued, based on mean field the- 
ory and projected wavefunction studies, that the U(l) 
uniform state which supports three spinon Fermi surfaces 
is an energetically viable candidate for the ground state 
of the s=l/2 hyper-kagome Heisenberg model. We have 
shown that it provides a reasonable overall description of 
the specific heat of Na^^Os over a broad temperature 
range T > 5K. Such spin liquids with spinon Fermi sur- 
faces have also been proposed recently for some quasi- 
two-dimensional frustrated magnets [IOJ. The spinon 
Fermi surfaces have direct physical implications for the 
low energy spin excitations. Specifically, from the mean 
field Hamiltonian, we can construct the gauge-invariant 
wavevectors connecting the surfaces, which dictate the 
wavevector dependence of triplet s = 1 states. 

At lower temperatures, there could be Fermi surface 
instabilities upon the inclusion of various small pertur- 
bations. For example, our analysis above did not include 
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FIG. 3: Comparison of the specific heat coefficient, C/T, of 
the U(l) uniform spin liquid state computed using RMFT 
(see text) with experimental data [lj on Na4lr30g after sub- 
traction of the phonon contribution. 



the effect of possible further neighbor interactions (small 
next neighbor antiferromagnetic exchange J' ~ 0.1— 0.2 J 
cannot be ruled out from our diagonalization study) . We 
have checked that extending our Zi ansatz to include 
next neighbor pairing terms, which would arise from fur- 
ther neighbor interactions in a mean field theory, leads to 
an unconventional pairing state with line nodes. These 
line nodes exist where the [110] plane (and symmetry re- 
lated planes) intersect the spinon Fermi surfaces. Such a 
line-node state would lead to a low temperature specific 
heat C ~ T 2 . As we mentioned, however, any Zi spin liq- 
uid state would most likely undergo a phase transition to 
the higher temperature paramagnetic phase, which was 
not observed in the experiment above 0.5-K\ Moreover, 
from Fig (31 it seems that one need not appeal to such 
a paired state to explain the linear-T behavior of C/T 
between 5K — 20K. If such a transition does exist, it is 
likely to occur at much lower temperature. 

Finally, it has been recently argued that the weak tem- 
perature dependence of \ a t low T and the specific heat 
C ~ T 2 cannot be reconciled unless spin orbit interac- 
tions are present 12j. It was shown that, despite rather 



strong atomic spin orbit coupling on Ir, the effective spin 
model is likely still of Heisenberg type with the dominant 
effect of spin-orbit induced Dzyaloshinskii-Moriya (DM) 
corrections. For sufficiently small DM (relative to J), 
our results for the energetics and specific heat would re- 
main unchanged. However, the DM coupling can strongly 
effect the spin susceptibility, especially in the line-node 
state. Such effects could potentially bring the suscep- 
tibility of the line-node state, which naively behaves as 
X(T) ~ T, into better agreement with experiment. The 
clarification of these issues is a promising direction for 
future research. 
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Note Added: During the final stages of the preparation 
of this manuscript, a paper dealing with related issues 
has appeared on arXiv [l3j. 
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